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ABSTRACT 

Magnetic fields in the early universe can significantly alter the thermal evolu- 
tion and the ionization history during the dark ages. This is reflected in the 21 cm 
line of atomic hydrogen, which is coupled to the gas temperature through colli- 
sions at high redshifts, and through the Wouthuysen-Field effect at low redshifts. 
We present a semi-analytic model for star formation and the build-up of a Lyman 
a background in the presence of magnetic fields, and calculate the evolution of 
the mean 21 cm brightness temperature and its frequency gradient as a function 
of redshift. We further discuss the evolution of linear fluctuations in temperature 
and ionization in the presence of magnetic fields and calculate the effect on the 
21 cm power spectrum. At high redshifts, the signal is increased compared to the 
non-magnetic case due to the additional heat input into the IGM from ambipolar 
diffusion and the decay of MHD turbulence. At lower redshifts, the formation of 
luminous objects and the build-up of a Lyman a background can be delayed by a 
redshift interval of 10 due to the strong increase of the filtering mass scale in the 
presence of magnetic fields. This tends to decrease the 21 cm signal compared 
to the zero-field case. In summary, we find that 21 cm observations may become 
a promising tool to constrain primordial magnetic fields. 

Subject headings: atomic processes — magnetic fields — turbulence — cosmology: 
theory — cosmology: early universe 
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1. Introduction 

Observations of the 21 cm fine struc- 
ture line of atomic hydrogen have the po- 
tential to become an important means of 
studying the universe at early times, dur- 
ing and even before the epoch of reion- 
ization. Th is possibility was sug gested 
originally by IPurcell fc Field! (119561 ). and 
significant process in instrumentation and 
the development of radio telescopes has 
brought us close to the first observa- 
tions from radio telescopes like LOFAR0. 
While one of the main purposes is to 
increase our unde rstanding of co s molog - 



ical reionization ( Ide Bruyn et al.l 120021 ) 



a number of further exciting applica- 
tions have been suggested in t he mean 
time. lLoeb &: Zaldarriagal (12004 ) demon- 
strated how 21 cm measurements can 
probe the thermal evolution of the IGM 
at a much earlier time, at redshi fts of 



200. iBarkana fc Loebl (120051 ) sug- 



gested a method that separates physical 
and astrophysical effects and thus allows 
to pro be the physics of the ear l y uni- 



verse. iFurlanetto. Oh fc Pierpaolil (120061 ) 
showed how the effect of dark matter an- 
nihilation and decay would be reflected in 
the 21 cm line, and effects of primordial 
magnetic fields have been considered by 
Tashiro fc Sugiyamal (120061 ). 



Indeed, primordial magnetic fields can 
affect the early universe in various ways. 
The thermal evolution is significantly al- 



Schleicher. Banerjee fc Klessenll2008l ). iKim. Olindo fc Rosn 
(119961 ) calculated the effect of the Lorentz 
force on structure formation and showed 
that additional power is present on small 
scales in the presence of primordial mag- 
netic fields. It was thus suggested that 

reionization occurs earlier in the presence 

of primordial magnetic field s (ISethi fc Subramanian 
20051 : iTashiro fe Sugivamal l2006h. How- " 



ever, as pointed out by iGnedinl (120001 ). the 
characteristic mass scale of star forming 
halos, the so-called filtering mass, increases 
signficantly when the temperature is in- 
creased. For comoving field strengths of ~ 
1 nG, we found that the filtering mass scale 
is shifted to scales where the power spec- 
trum is ess entially independent of the mag- 
netic field ( jSchleicher. Banerjee fc Klessen 
20081 ). Reionization is thus delayed in 



the presence of primordial magnetic fields. 
We further found upper limits of the or- 
der 1 nG, based on the Thomson scat- 
te ring optical depth measured by WMAP 
5 (iKomatsu et al.ll2008l : iNolta et al.l 120081 ) 
and the r equirement that rei onization ends 
at z ~ 6 dBecker et al.ll200lh . 

The presence of primordial magnetic 
fields can have interesting implications on 
first star formation as well. In the absence 



of ma gnetic fields. lAbel. Bryan &: Norman 
(120021 ) and iBromm fc Larson! (12004 ) sug- 
gested that the first stars should be very 
massive, perhaps with ~ 100 solar masses. 



Clark. Glover fc Klessen! (120081 ) and lOmukai. Schneider fc F. 



decaying MHD turbulence ( 


Sethi & Subramanian 


2005 




Sethi, Nath & Subramanian 


2008 





^lttpi/Zwww. lofar.org 



(120081 ) argued that in more massive and 
perhaps metal-enriched galaxies, frag- 
mentation should be more effective and 
lead to the formation of rather low-mass 
st ars , due to a stage of efficien t coo l- 
ing (jOmukai. Schneider &: Haimanl 120081 ). 
Magnetic fields may change this picture 
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and reduce the stell ar mass by triggerin g 
jets and outflows flSilk fc Langerl 120 06). 
However, simulations by (IMachida et al. 



20061 ) show that the change in mass is of 
the order 10%. 



21 cm measurements can try to adress 
primordial magnetic fields in two ways: 
During the dark ages of the universe, at 
redshifts z ~ 200 well before the forma- 
tion of the first stars, the spin temperature 
of hydrogen is coupled to the gas temper- 



drogen atoms ( 


Allison & Dalearno 


1969; 


Zveelman 


2005T) and free electrons ( 


Smith 


1966: 


Furlanetto & Furlanetto 


2007), con- 



stituting a probe at very early times. 
While collisional de-excitation becomes in- 
efficient due to the expansion of the uni- 
verse, the first stars will build up a Lyman 
a background that will cause a deviation 
of the spin temperature from the radiation 
temp erature by the Wouthuysen-Fie ld ef- 
fect flWouthuvsenl Il952l : iFieldl ll95Sh . As 



primordial magnetic fields may shift the 
onset of reionization, the onset of this cou- 
pling constitutes an important probe on 
the presence of such fields. We adress 
these possibilities in the following way: In 
§0 we review our treatment of the IGM in 
the presence of primordial magnetic fields. 
The evolution of the 21 cm background and 
the role of Lyman a photons is discussed in 
§|3j The evolution of linear perturbations 
in temperature and ionization is calculated 
in HI The results for the power spectrum 
are given in §51 and the results are further 
discussed in 361 



2. The evolution of the IGM 

As indicated in the introduction, pri- 
mordial magnetic fields can have a strong 
impact on the evolution of the IGM be- 
fore reionization. Once the first luminous 
objects form, their feedback must also be 
taken into account. In this section, we 
review the basic ingredients of our treat- 
ment of the IGM between recombination 
and reionization. We refer the interested 



reader to ISchleicher. Banerjee fc Klessen 
( 120081 ) for more details. 



2.1. The RECFAST code 

Our calculation is based on a modified 



version of t he RECFAST codcl ( ISeager. Sasselov fc Scott 
20001 . Il999l ) that calculates recombination 
and the freeze-out of electrons. The cal- 
culation of helium recombination was re- 



cently updated by IWong. Moss fc Scott 
( 120071 ) . We have extended this code 
by including a model for reionization in 
the context of primordial magnetic fields 
Schleicher. Banerjee &: Klessenl ( 120081 ) . The 
equation for the temperature evolution is 
given as 



dT 

dz 



8a T a R Tr 



<T-T r 



3H(z)(l + z)m e c 1 + /hc + x, 
IT 2(Lh cat — -^cooi) 
l + z ~ 3nk B H(z)(l + z)' 

where Lh ea t is the new heating term (see 
§221 and §GL2D, L cool the cooling by 
Lyman a emission, bremsstrahlung and re- 
combinations, ot is the Thomson scatter- 
ing cross section, a R the Stefan-Boltzmann 
radiation constant, m e the electron mass, 
c the speed of light, k R Boltzmann's con- 
stant, n the total number density, x e = 



(1) 



2 http: / /www. astro. ubc.ca/ people / scott /recfast .html 
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Fig. 1. — Evolution of gas temperature 
(upper panel) and ionized fraction (lower 
panel) in the medium that was not yet af- 
fected by UV feedback. The details of the 
models are given in Table [TJ 
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Fig. 2. — Evolution of the filtering mass 
scale for different comoving magnetic field 
strengths. 



ra e /?T,H the electron fraction per hydrogen 
atom, H(z) is the Hubble factor and /n e is 
the number ratio of He and H nuclei, which 
can be obtained as /n e = Y p /A{l — Y p ) from 
the mass fraction Y p of He with respect to 
the total baryonic mass. Eq. p]) describes 
the change of temperature with redshift 
due to Compton scattering of CMB pho- 
tons, expansion of the universe and addi- 
tional heating and cooling terms that are 
described in the following subsections in 
more detail. The evolution of the ionized 
fraction of hydrogen, x p , is given as 



dz 



[x e x p n n a n - Pn(l 



x p )e 



X 



H(z)(l +z)[l + K U (A R + (3 H )n u (l - z P )] 
[1 + K H A H n H (l - x P )] 



H(z)(l + z) 



fk 



(2) 



where the ionization fraction is deter- 
mined by radiative recombination, photo- 
ionization of excited hydrogen atoms by 
CMB photons in an expanding universe 
and collisional ionization, as well as X-ray 
feedback. Here, nn is the number density 
of hydrogen atoms and ions, h p Planck's 
constant, f lon describes ionization from X- 
rays and UV photons (see §2.31) . and the 
parametrized case B recombination coeffi- 
cient for atomic hydrogen an is given by 



a H = F x 10 13 at , cm 3 s 
1 + ct d 



(3) 



with a = 
0.6703, d 



4.309, b = 
-- 0.5300 and 



-0.6166, c 
t = T/10 4 



K 



(jPeauignot et al.lll99ll : iHummer fc Storey 
. This coefficient takes into account 
that direct recombination into the ground 
state does not lead to a net increase of 
neutral hydrogen atoms, since the photon 
emitted in the recombination process can 
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ionize other hydrogen atoms in the neigh- 
bourhood. 

The fudge factor F = 1.14 serves to 
speed up recombination and is determined 
from comparison with the multilevel-code. 
The photoionization coefficient /% is cal- 
culated from detaile d balance at high red- 



shifts a s desc ribed by lSeager. Sasselov & Scott 
(|2000L Il999h . Once the ionized fraction 
drops below 98%, it is instead calcu- 
lated from the photoionization cr o ss sec - 
tion given by ISasaki fc Takaharal (119931 ). 
as detailed balance may considerably over- 
estimate the photoionization coefficient in 
case of additional energy input at low red- 
shift. The wavelength Xu,2 P corresponds 
to the Lyman-a transition from the 2p 
state to the Is state of the hydrogen atom. 
The frequency for the two-photon transi- 
tion between the states 2s and Is is close 
to Lyman-a and is thus approximated 
by ^h,2s = c/Ah,2 P (he., the same aver- 
aged wavelength is used). Finally, Ah = 
8.22458 s _1 is the two-photon rate for 
the tr ansition 2s-ls according to iGoldman 
fll989f ) and K H = \^/[S7rH{z)] the cos- 
mological redshifting of Lyman a pho- 
tons. To take into account UV feedback, 
which can not be treated SIS db background 
but creates locally ionized bubbles, we 
further calculate the volume-filling factor 



medium to temperatures above 10 4 K be- 
fore it was ionized. 

2.2. Heating due to primordial mag- 
netic fields 

The presence of magnetic fields leads 
to two different contributions to the heat- 
ing rate, one coming from ambipolar dif- 
fusion and one resulting from the de- 
cay of MHD turbulence. In the first 
the contribution can be calculated 



as ( 


Sethi & Subramanian 




2005; 


Cowline 


1956 




Schleicher. Baneriee 


& Klessen 


2008) 



7 ambi 



Pn 



B [ 



16TT 2 -fp 2 b p l L 2 



(4) 



Here, p n , pi and pb are the mass densi- 
ties of neutral hydrogen, ionized hydrogen 
and all baryons. The ion-neutral coupling 
coefficient is calculated using the updated 
zero drift ve locity moment u m tra nsfer co- 
efficients of IPinto fc Gailil fj2008h for col- 
lisions of H + with H and He. The coher- 
ence length L is estimated as the inverse of 



by 


— . , n . , „ . n , _ TT|| 

Jedamzik. Katalinic & Olinto 




1998 




Subramanian k, Barrowlll998; 


Seshadri & Subramanian 


2001 


1 as 



234 Mpc" 



B, 







prescriptions of ( 


Shapiro & Giroux 1987 




Haiman & Loeb 


1997 


Barkana & Loet 


2001;lLoeb & Barkana 


2001: Choudhurv & Ferrara 


2005: Schneider et al. 


2006). We as- 


which 



X 



Q b h 2 



1/2 



lO" 9 G 
^ 1/4 



f2 m 
(L3 



0.7 



1/4 



(5) 



sume that the temperature in the ionized 
medium is the maximum between 10 4 K 
and the temperature in the overall neu- 
tral gas. The latter is only relevant if pri- 
mordial magnetic fields have heated the 



the magnetic field are damped out during 
recombination. The comoving magnetic 
field is denoted as B = B/(l + z) 2 . 

For decaying MHD turbulence, we adopt 
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the p rescription of ISethi fc Subramanian 
f|2005h . 



^decay 



B (t) 2 3a [ln(l + t d /U)fH(t) 



8tt 2 [ln(l + U/U) + ln(t/ti)] 5+1 

(6) 

where i is the cosmological time at red- 
shift z, td is the dynamical timescale, U 
the time where decay starts, i. e. after the 
recombination epoch when velocity per- 
turbations are no longer damped by the 
large radiative viscosity, Zi is the corre- 
sponding redshift. For a power spectrum 
of the magnetic field with power-law index 
a, implying that the magnetic field scales 
with the wavevector k to the power 3 + a, 



3)/(a + 5) ( 


Olesenl Il997: ShiromizJ 1998: 


Christensson. Hindmarsh & Brandenburd 


2001; 


Banerjee & Jedamzik 


2003). In the 



generic case, we expect the power spectrum 
of the magnetic field to have a maximum at 
the scale of the coherence length, and the 
heat input by MHD decay should be deter- 



ing to larger scales ( 


Miiller & Biskamp 


2000; 


Christensson. Hindmarsh & Brandenburg 


2001: 


Banerjee & Jedamzik 2003 


■ 


2004 





We thus adopt a = 3 for the calculation. 
We estimate the dynamical timescale as 
td = L/va, where va = B / is the 
Alven velocity and p\> the baryon mass 
density. The evolution of the magnetic 
field as a function of redshift can be de- 
termined from the magnetic field energy 
density E n = B 2 /87i : which evolves as 



(ISethi fe Subramanian! 120051 ) 



dE 



B 



dt 



—4H(t)EB — -^ambi ~ ^deeay (7) 



For the models given in Table [U the evolu- 
tion of temperature and ionization in the 



overall neutral medium, i. e. the gas that 
was not yet affected by UV photons from 
reionization, is given in Fig. [TJ 

For comoving fields weaker than 0.02 nG, 
the gas temperature is closely coupled to 
the CMB via Compton-scattering at red- 
shifts z > 200. At lower redshifts, this 
coupling becomes inefficient and the gas 
cools adiabatically during expansion, un- 
til heating due to X-ray feedback becomes 
important near redshift 20. The ionized 
fraction drops rapidly from fully ionized at 
z ~ 1100 to ~ 2 x 10 _4 , until it increases 
again at low redshift due to X-ray feed- 
back. For larger magnetic fields, the addi- 
tional heat input becomes significant and 
allows the gas to decouple earlier from the 
CMB. For comoving fields of order 0.5 nG 
or more, the gas reaches a temperature 
plateau near 10 4 K. At this temperature 
scale, collisional ionizations become im- 
portant and the ionized fraction in the gas 
increases, such that ambipolar diffusion 
heating becomes inefficient. The onset of 
X-ray feedback is delayed to the higher fil- 
tering mass in the presence of magnetic 
fields, which is discussed in the next sec- 
tion. 

2.3. The filtering mass scale and 
stellar feedback 

The universe becomes reionized due to 
stellar feedback. We assume here that the 
star formation rate (SFR) is proportional 
to the change in the collapsed fraction / co u, 
i. e. the fraction of mass in halos more mas- 
sive than m m - m . This fraction is g i ven b y 



the formalism of Press & Schechtcr ( 



In ISchleicher. Banerjee fc Klessenl ( 12008 



1974). 



we have introduced the generalized filter- 
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mg mass m F B} given as 
3 



M 



2/3 
F.B 



(l 



da'M 2 /\a! 



(r 






5 



where a = 1/(1 + z) is the scale factor and 
M 9 is the maximum of the thermal Jeans 
mass Mj and the magnetic Jeans mass 
Mf , the mass scale below which magnetic 



tional colapse ( 


Sethi & Subramanian 


2005; 


Subramanian & Barrow 


1998). The con- 



cept of the filtering mass, i. e. the halo 
mass for which the baryonic and dark 
m atter evolution c an de couple, go es back 
to iGnedin fc Huil (119981 ) , see also iGnedin 



fl2000h . To take into account the back re- 



action of the photo-heated gas on struc- 
ture formation, the Jeans mass is cal- 
culated from the effective temperature 
T e ff = QunTmax + (1 - Qhtt)T, where 
T mnx = max(10 4 K,T) (ISchneider et al. 



20061 ). The evolution of the filtering mass 
for different magnetic field strengths is 
given in Fig. 12.11 

When the comoving field is weaker than 
0.03 nG, the filtering mass is of order 
10 5 Mq at high redshift and increases up 
to 10 7 Mq during reionization, as UV and 
X-ray feedback heats the gas temperature 
and prevents the collapse of baryonic gas 
on small scales. Stronger magnetic fields 
provide additional heating via ambipolar 
diffusion and decaying MHD turbulence, 
such that the filtering mass is higher from 
the beginning. This implies that the first 
luminous objects form later, as the gas in 
halos below the filtering mass cannot col- 
lapse. This effect is even stronger when 
the comoving field is stronger than 0.3 nG, 
as the magnetic Jeans mass then domi- 
nates over the thermal Jeans mass and 



increases the filtering mass scale by fur- 
ther orders of magnitude. As discussed in 



Schleicher. Banerjee &: Klessenl (120081 ) , the 



feedback on structure formation via the fil- 
tering mass can delay reionization signifi- 
cantly, which in turn allows to calculate 
upper limits on the magnetic field strength 
due to the measured reionization optical 
depth. 

Consequently, we define the lower limit 
to halo masses that can form stars as 
m min = max(M FjB , 10 5 M Q ), where 10 5 M© 
is the minimal mass scale for which baryons 
c an cool efficiently, as found in simulations 
of lGreif et al.l (120081 ). We take into account 
X-ray feedback that can penetrate into the 
IGM due to its long mean-free path, as well 
as feedback from UV photons produced in 
the star forming regions. As w e showed 



m iSchleicher. Banerjee fc Klessenl ( 120081 ). 
Population II stars can not reionize the 
universe efficiently. It is thus reasonable 
to assume that the first luminous sources 
were indeed massive Pop. Ill stars. While 
this may no longer be true once chemical 
and radiative feedback becomes efficient 
and leads to a different mode of star forma- 
tion, such a transition seems not important 
for the main purpose of this work, which is 
to demonstrate how primordial magnetic 
fields shift the onset of reionization and 
the epoch where a Lyman a background 
bu ilds up initially. We thus adopt model B 
of ISchleicher. Banerjee fc Klessenl ( 120081 ). 
which uses an escape fraction of 100% if the 
virial temperature is below 10 4 K, and 10% 
for the atomic cooling halos with higher 
virial temperature. It further assumes that 
4 x 10 4 UV photons are emitted per stellar 
baryon during the lifetime of the star. This 
choice is in particular motivated by simula- 
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tions of IWhalen. Abel fc Norman! (12004 ), 
that show that minihalos are easily photo- 
evaporated by UV feedback from Pop. 
Ill stars, yielding escape fr actions close 
to 100%. A new study by IWise fc Cen 



(120081 ) shows that the escape fraction may 
be larger than 25% even for atomic cool- 
ing halos. However, their study assumes 
purely primordial gas, which is unlikely in 
such systems, and we expect that studies 
which take into account metal enrichment 
will find lower escape fractions, perhaps of 
the order 10% as suggested here. 



3. The 21 cm background 

The 21 cm brightness temperature de- 
pends on the hyperfinestructure level pop- 
ulations of neutral hydrogen, which is de- 
scribed by the spin temperature T spin . In 
this section, we review the physical pro- 
cesses that determine the spin tempera- 
ture, discuss the build-up of a Lyman-a 
background that can couple the spin tem- 
perature to the gas temperature at low red- 
shift, and discuss various sources of 21 cm 
brightness fluctuations. 

3.1. The spin temperature 

The observed 21 cm brightness temper- 
ature fluctua tion can then be conveniently 
expressed as (IFurlanetto. Oh fc Briggsll2006l ) 



ST h 



27x H (l + 5) 



Tft-Tr 



Q h h 2 



0.15 l + z 



x 



0.023/ \a m h 2 10 



H(z)/(l + z) 



dv\\/dr\ 



where 5 is the fractional overdensity, 
ih = 1 — ^eff the effective neutral fraction, 
x e f{ = Qh+ + (1 — Qu+) x p the effective ion- 
ized fraction, and Q m the cosmological 



1/2 



mK, (9) 
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Fig. 3. — CMB, gas and spin temperature 
in the gas unaffected by reionization, for 
a star formation efficiency /* = 10~ 3 , in 
the zero field case (upper panel) and for 
a field strength of 0.8 nG (lower panel). 
In the absence of magnetic fields, gas and 
CMB are strongly coupled due to efficient 
Compton scattering of CMB photons. At 
about redshift 200, the gas temperature de- 
couples and evolves adiabatically, until it 
is reheated due to X-ray feedback. The 
spin temperature follows the gas tempera- 
ture until redshift ~ 100. Collisional cou- 
pling then becomes less effective, but is re- 
placed by the Wouthuysen-Field coupling 
at z ~ 25. In the presence of magnetic 
fields, additional heat goes into the gas, it 
thus decouples earlier and its temperature 
rises above the CMB. Wouthuysen-Field 
coupling becomes effective only at z ~ 15, 
as the high magnetic Jeans mass delays the 
formation of luminous sources significantly. 



3 1000 r 




Wavelength [in] 
10 




Wavelength [ml 
10 



100 




Fig. 4. — The mean evolution in the over- 
all neutral gas, calculated for the models 
described in Table [TJ Top: The spin tem- 
perature. Middle: The mean brightness 
temperature fluctuation. Bottom: The fre- 
quency gradient of the mean brightness 
temperature fluctuation. 



density parameters for baryonic and to- 
tal mass, and dvu/dru the gradient of the 
proper velocity along the line of sight, in- 
cluding the Hubble expansion. The spin 
temperature T spin is given as. 



spin 



(10) 



where T c is the colour temperature, x c 



x" + x e c the collisional coupling co- 
efficient, given as the sum of the cou- 
pling coefficient with respect to hydrogen 
atoms, xf H = Ki_ nnT^/ A w T r , and the 
coupling coefficient with respect to elec- 
trons, x e c H = Kl_ n e T*/A w T r . Further, 
Uh denotes the number density of hydro- 
gen atoms, T* = 0.068 K the temperature 
corresponding to the 21 cm transition and 
A w = 2.85 x 10~ 15 s" 1 the Einstein co- 
efficient for radiative de-excitation. For 
the collisional de-excitation rate through 
hydrogen, K?_ n , w e use the recent calcu- 



lation of IZvgelmanl f|2005l), and we use 
the new rate of Furlanetto fc Furlanetto 



( 120071 ) for the collisional de-excitation rate 
through electrons, Ki_ Q . When both the 
spin and gas temperature are significantly 
larger than the spin exchange temperature 
T se = 0.4 K, as it is usually the case, the 
colour temperature is given as 



T r = T 



1 + T se /T 
1 T se /T S pi n 



:n) 



Spin and colour temperature thus depend 
on each other and need to be solved for si- 
multaneously. The Wouthuysen-Field cou- 
pling coefficient is given as 



"c 
10/ 



(12) 



where J c = 1.165x 10- 10 (l+z)/20 cm" 2 s" 1 Hz 
and J a is the proper Lyman a photon in- 



sr 



9 



tensity, given as 



"max 



JJz) 



n=2 



n 



2max , ,(1 + z) 2 
02 - 



4tt H(z' 
(13) 

(see also iBarkana fc Loebl 120051 ). / rec is 
the r ecycling probability of Lyman n pho- 
tons (IPritchard fc Furlanettdl2006l ; iHirata 
20061 ). H(z) the Hubble function at red- 
shift z, c the speed of light, n max the high- 
est considered Lyman n resonance, z max 
the highest redshift at which star forma- 
tion might occur and e(v, z) the comov- 
ing emissivity. We adopt n max = 30 and 



40. The suppresion factor S 



be calculated as ( Chuzh ov fc Shapiro 



Furlanetto fc Pritachardl 120061 ) 



can 



2006 



S a = exp (-0.803T~ 2/3 (lO- 6 r GP ) 1/3 

(14) 

from the optical depth for Lyman a pho- 
tons, 



r G p ~ 3 x 10 x hi 



1 + z 
7 



3/2 



(15) 



Examples for the evolution of the gas, spin 
and CMB temperature in the zero-field 
case and for a magnetic field of 0.8 nG are 
given in Fig. [3J 

3.2. The build-up of a Lyman-a 
background 

In the dark ages, the universe is trans- 
parent to photons between the Lyman 
a line and the Lyman limit. Luminous 
sources can thus create a radiation back- 
ground between these frequencies that is 
redshifted during the expansion of the uni- 
verse. When a photon is shifted into a 
Lyman line, it is scattered and can couple 



-e f„ 



the spin temperature to the gas temper- 
ature as des cribe d above. As shown b y 



Hjratal (120061 ) and lBarkana fc Loebl (120051 ) 



a population of Pop. Ill stars can already 
produce a significant amount of Lyman 
a radiation around redshift 25. We as- 
sume that emission is proportional to the 
star formation rate, which is assumed to 
be proportional to the change in the col- 
lapsed fraction / co n. It thus depends on the 
evolution of the generalized filtering mass 
in the pre sence of magnetic fields. As dis- 
cuss ed in iFurlanetto. Oh fc Briggsl (120061 ) 
and IBarkana &: Loebl (120051 ) , the comoving 
emissivity is given as 



e z/, z) 



f*n b e b (v 



df, 



coll 



dt 



(16) 



where n b is the comoving baryonic number 
density and e b {y) the number of photons 
produces in the frequency interval v±dv/2 
per baryon incorporated in stars. The 
emissivity is approximated as a power law 
€b(v) oc z/ Qs_1 , and the parameters must 
be chosen according to the model for the 
stellar population. As discussed above, 
we assume that the first stars are massive 
Pop. HI stars, and a dopt the model of 
Barkana k Loebl (|2005h with N a = 2670 



photons per stellar baryon between the Ly- 
man a and the Lyman f3 line and a spectral 
index a s = 1.29 

In addition, secondary effects from X- 
ray photons may lead to excitation of 
hydrogen atoms and the production of 
further Lyman a photons. We calcu- 
late the production of Lyman a photons 
from the X-ray background for consis- 



2006 



given as ( 


Chen & Miralda-Escude 


Chuzhov. Alvarez & Shapiro 




2006; 
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Furlanetto. Oh k BriggsH2006T ) 

A X-ray n HZQ f ( fx,caW f* coll /dz 

A Xa = omsM [^-qj-q^- 

x(^) 3 - 

We adopt fx = 0.5, a nd calculate /y.roii 
from t he fit-formulae of IShull fc Steenberg 
(119851) . (Notation: Please recall that / co n 
denotes the fraction of collapsed dark mat- 
ter, while /x.coii denotes the fraction of 
absorbed X-ray energy which goes into 
X-rays. The parameter fx provides an 
overall normalization for the correlation 
between the star formation rate and the 



X-ray luminosity (IKoyama et al.l Il995l ).) 
For the models described in Table [TJ we 
have calculated the evolution of the spin 
temperature and the mean 21 cm bright- 
ness fluctuation. From an observational 
point of view, it might be more reason- 
able to focus on the frequ ency gradient of 



the m ean 21 brightness ([Schneider et al. 



20081 ) . As the frequency dependence of the 
foreground is known, such an analysis may 
help to distinguish between the foreground 
and the actual signal. Results for the spin 
temperature, the mean 21 cm brightness 
fluctuation and its frequency gradient are 
shown in Fig. HI 

At redshifts above 200, the spin tem- 
perature is coupled closely to the gas tem- 
perature due to collisions with hydrogen 
atoms. At lower redshifts, the coupling 
becomes inefficient and the spin temper- 
ature evolves towards the CMB tempera- 
ture, until structure formation sets in and 
the spin temperature is again coupled to 
the gas temperature via the Wouthuysen- 
Field effect. The point where this happens 
depends on the magnetic field strength, 



which influences the filtering mass and 
thus the mass scale of halos in which stars 
can form. Especially in the presence of 
strong magnetic fields, this delay in struc- 
ture formation also makes the production 
of Lyman a photons less efficient and the 
departure from the CMB temperature is 
only weak. 

This bevahiour is reflected in the evo- 
lution of the mean brightness temperature 
fluctuation. For weak magnetic fields, the 
gas temperature is colder than the CMB 
because of adiabatic expansion, so the 21 
cm signal appears in absorption. At red- 
shift 40 where the spin temperature ap- 
proaches the CMB temperature, a maxi- 
mum appears as absorption goes towards 
zero. Comoving fields of 0.05 nG add only 
little heat and essentially bring the gas 
and spin temperatures closer to the CMB, 
making the mean brightness temperature 
fluctuation smaller. For stronger fields, 
the 21 cm signal finally appears in emis- 
sion. Due to very effective coupling and 
strong temperature differences at redshifts 
beyond 200, a pronounced peak appears 
there in emission. At redshifts between 
20 and 10 which are more accessible to 
the next 21 cm telescopes, magnetic fields 
generally reduce the expected mean fluc- 
tuation, essentially due to the delay in the 
build-up of a Lyman a background. 

The evolution of the frequency gradient 
in the mean brightness temperature fluctu- 
ation essentially shows strong minima and 
maxima where the mean brightness tem- 
perature fluctuation changes most signif- 
icantly. In the case of weak fields, this 
peak is near redshift 30 when coupling via 
the Wouthuysen-Field effect becomes ef- 
ficient. For stronger fields, this peak is 
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shifted towards lower redshifts, thus pro- 
viding a clear indication regarding the de- 
lay of reionization. 

Scattering of Lyman a photons in the 
IGM does not only couple the spin tem- 
perature to the gas temperature, but it 
is also a potential source of heat. Its 
effect on the mean temperature evolu- 



authors ( 


Madau. Meiksin & Rees 


1997 




Chen & Miralda-Escudel 


20061: Chen&Miralda-Escude! 


2004 


Chuzhov & ShaDiroll2006l: 


Ciardi, Bianchi & Ferrara 


2002). While a detailed treatment requires 



Model B Q [nG] /* 



to separately follow the mean tempera- 
ture of hydrogen and deuterium, a good 
estimate can be obtained using t he a 
proach of iFurlanetto. Oh k, Briggj (120061 
Furlanetto fc Pritachardl fl2006h . i. e. 



2 L 



heat,Lya 



0.8 x a ( 10 



3n n k B TH{z) T 4 / 3 S a \1 + z 



(18) 



For star formation efficiencies of the order 
0.1%, as they are adopted throughout most 
of this work, Lyman a heating leads to a 
negligible temperature change of the order 
1 K. 

3.3. 21 cm fluctuations 

Fluctuations in the 21 cm brightness 
temperature can be caused by the rela- 
tive density fluctuations 5, relative tem- 
perature fluctuations 5t and relative fluc- 
tuations in the neutral fraction 5h- Fluc- 
tuations in the Lyman a background can 
be an additional source of 21 cm fluctu- 
ations (jBarkana fc Loebl 120051 ). but a de- 
tailed treatment of these fluctuations is be- 
yond the scope of this work. These fluc- 
tuations are only important in the early 
stage of the build-up of such a background. 
Once x a is significantly larger than unity, 



1 





0.1% 


2 


0.02 


0.1% 


3 


0.05 


0.1% 


4 


0.2 


0.1% 


5 


0.5 


0.1% 


6 


0.8 


0.1% 


7 


0.8 


1% 



Table 1: A list of models for different co- 
moving magnetic fields and star formation 
efficiencies, which are used in several fig- 
ures for illustrational purposes. We give 
the comoving magnetic field Bq and the 
star formation efficiency /*. For illustra- 
tion purposes, all models assume a pop- 
ulation of massive Pop. Ill stars. The 
amount of Lyman a photons produced per 
stellar baryon would be larger by roughly 
a factor of 2 if we were to assume Pop. 
II stars. Assuming that the same amount 
of mass goes into stars, the coupling via 
the Wouthuysen-Field effect would start 
slightly earlier, but the delay due to mag- 
netic fields is still more significant. 
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the spin temperature is coupled closely 
to the gas temperature, and small fluctu- 
ations in the background will not affect 
the coupling. We adopt t he treatment of 
Loeb fc Zaldarriagal (12004 ) to calculate the 
fractional 21 cm brightness temperature 
perturbation, given as 

hi{k) = ((3 + p 2 )5 + (3 H 5 H + (3 T 5 T , (19) 

where p is the cosine of the angle be- 
tween the wavevector k and the line 
of sight, with the expansion coefficients 



scale dependence ( 


Barkana & Loeb 


2005; 


Naoz & Barkana 


2005). However, on the 



large scales which might ultimately be ob- 
serveab le, the growing density mode dom- 



inates (IBharadwaj fc Alill2004l ). The tern 



perature and ionization fluctuations 5t and 
5i are thus related to the density fluctua- 
tions 5 by 5t = gr(z)5 and 5i = gi(z)5 
with redshift-dependent coupling factors 
gi(z) and gr{z)- The time evolution of 
qr and o,- is well-known in the absence 



( Ba rkana fc Loebll2005l : iFurlanetto. Oh fc Brk 
2006h 



of magnetic fields 


(Naoz & Barkana 


2005; 


Bharadwaj & Ali \ 


>004 


) and has also been 



studied for the case of dark matter annihi- 



1 + 



^tot (1 + X t 



(20) 



1 + 



x 



HH 



X 



at J 
eH 



^tot (1 + 



(21) 



dlnx, 



T — T r x tot (l + x tot ) din 



The coefficient x to t is given as the sum of 
the collisional coupling coefficient x c and 
the Wouthuysen-Field coupling coefficient 
x a . 

4. The evolution of linear perturba- 
tions in the dark ages 

To determine the fluctuations in the 
21 cm line, we need to calculate the evo- 
lution of linear perturbations during the 
dark ages. We introduce the relative per- 
turbation of the magnetic field, 5b, and 
the relative perturbation to the ionized 
fraction, <5j, which is related to the rela- 
tive perturbation in the neutral fraction, 
<5h, by Si = -S n (l - Xi)/xi. In the gen- 
eral case, different density modes trigger 
independent temperature and ionization 
fluctuations. This introduces a non-trivial 



lation and decay (IFurlanetto. Oh fc Pierpaoli 
20061 ). We further introduce the relative 
fluctuation in the magnetic field, 8b = gB<5- 
To extend the previous analysis to the sit- 
uation considered here, we must consider 
the effect of the new terms due to ambipo- 
lar diffusion heating and decaying MHD 
turbulence in eq. ([T|). As it has a compli- 
cated dependence on different quantities, 
we will keep this discussion rather generic, 
so that it can also be applied to other sit- 
uations as well. 

Let us consider a source term of the 
form 



5T 



f{p,T,x h B). (23) 



Here, f(p,T,Xi,B) corresponds to the 
term describing ambipolar diffusion heat- 
ing and decaying MHD turbulence in Eq. 
(CE]). When the quantities p, T, Xi and B 
are perturbed, it is straightforward to show 
that this introduces further terms for the 
evolution of the temperature perturbation, 
which are of the form 



55t 

77 



_! df df df 

i5p dp- + T5T df + 5 ^dx- t 
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+ 5 B B 



df_ 
dB 



(24) 



These new terms describe the change in 
the relative temperature fluctuation, de- 
pending on the additional source term 
f(p,T,Xi,B) and its derivatives with re- 
spect to gas density, temperature, ionized 
fraction and magnetic field. The behaviour 
is dominated by the properties of the am- 
bipolar diffusion heating term. The time 
evolution of the coupling factors gr and gi 
is thus given as 

dgr _ ^t-2/3 + Xi g T T r - gi(T r - T) 



dz 



dgi 
dz 



1 + z 77it 7 T(l + z)H(z) 
1( df df df 

+ {p d-p +T9T df + 9 ^dx- 

+ gnB^-fgr), (25) 

9i | oiHXinH (1 + 9i + a- H 9: 



1 + z 



[l+z)H(z) 



The first term in Eq. ([25]) describes the 
evolution towards an adiabatic state in the 
absence of heating or cooling mechanisms. 
The second term describes the interaction 
with the CMB via Compton scattering, 
which in general drives the gas towards 
an isothermal state, as the heat input per 
baryon is constant. The further terms are 
those derived in Eq. (JMD- In Eq. Q2SD, the 
first term holds the relative fluctuation 8i 
constant in the absence of recombinations, 
while the second term describes the effect 
of hydrogen recombination, which tends to 
drive the corresponding coupling factor 
towards —1. Fluctuations in density and 
ionization are then anticorrelated. 

This system must be closed with an ad- 
ditional assumption for the fluctuation in 
the magnetic field. The dominant mecha- 
nism for energy loss is ambipolar diffusion, 



which scales with p^ 2 . It seems thus un- 
likely that the magnetic field will be dissi- 
pated in regions of enhanced density. We 
rather assume that the magnetic flux is 
approximately conserved. In the case of 
dynamically weak field strength, the cor- 
responding surface scales as p -2 ^ 3 , so the 
magnetic field scales as p 2//3 , thus gs ~ 
2/3. Ambipolar diffusion, the dominant 
new heating term, scales as p~ 2 . To re- 
late the change in energy to the change 
in temperature, it is further multiplied by 
3/(2nkB), such that in this case / oc p -3 . 
This behavior is not fully compensated by 
the other terms, thus heating is less effec- 
tive in overdense regions, and gr tends to 
drop below the adiabatic value of 2/3. For 
the models given in Table [lj the evolution 
of gx and ^ is given in Fig. [5J 

In the zero or weak magnetic field 
case, the coupling factor evolves from 
zero (efficient coupling to the CMB) to- 
wards 2/3 (adiabatic), until X-ray feed- 
back sets in. For comoving fields larger 
than 0.05 nG, there is an early phase 
around z ~ 500 where it evolves towards 
adiabatic behaviour, as the gas recombines 
more efficiently in these early times and 
higher densities, where the factor gi is still 
of order — 1. Thus, ambipolar diffusion in 
density enhancements is more efficient and 
effectively heats the gas at these redshifts. 
Below redshift 500, the relative fluctua- 
tion in ionization evolves towards zero, as 
the gas becomes thinner and hotter. That 
leads to a phase in which the relative tem- 
perature fluctuations evolves towards an 
isothermal state (gr = 0). At low red- 
shifts, X-ray feedback leads to a decrease 
of the coupling factor g?- Its behaviour 
is no longer adiabatic because of this ad- 
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ditional heat input. The increase in the 
mean ionized fraction due to secondary 
ionizations leads to a more negative gi. 

5. The 21 cm power spectrum 

We can now calculate the 21 cm power 
spectrum in the presence of primordial 
magnetic fields from eq. (TL9l . leading to 

P 21 (k,n) = 5T 2 ((3' + fi 2 ) 2 P ss (k), (27) 

where Pss{k) is the baryonic power spec- 
trum and 

P = P + f3 T g T - P H x i9i /(l - Xi ). (28) 

For simplicity, we average over the \x de- 
pendence and define P$s = Po/(l + z) 2 , 
such that 



Wavelength [m] 



100 



P2i{k) 



5T 2 



& 



12 



3 H 5 



Po(k) 



:i+^) 2 

(29) 

On the large scales considered here, the 
power spectrum of the 21 cm line is pro- 
portional to the baryonic power spectrum 
Po(k) f° r a given redshift, but the pro- 
portionality constant will evolve with red- 
shift and is independent of the wavenum- 
ber in the linear approximation. The ratio 
P21/P0 is shown in Fig. [5j It essentially re- 
flects the behaviour that was found earlier 
for the mean brightness temperature fluc- 
tuation. At high redshift, the amplitude 
of the 21 cm power spectrum is larger in 
the presence of additional heat from am- 
bipolar diffusion and decaying MHD tur- 
bulence. At low redshifts, it is smaller 
in the presence of magnetic fields, as the 
high filtering mass delays the build-up of 
a Lyman a background and leads only to 
a weak coupling of the spin temperature 
to the gas temperature. In addition, the 
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Fig. 5. — The evolution of the 21 cm 
large-scale power spectrum and the growth 
of large-scale fluctuations in temperature 
and ionization for the models given in Ta- 
ble [TJ Top: The evolution of the 21 cm 
power spectrum, normalized with respect 
to the baryonic power spectrum at red- 
shift 0. Middle: The evolution of the ratio 
qt = St/8- Bottom: The evolution of the 
ratio gi = 5 X J5. 
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point where this coupling starts is shifted 
towards lower redshifts. This provides two 
signatures that can provide strong indica- 
tions for the presence of primordial mag- 
netic fields. 

At high redshifts z > 50, the additional 
heat increases the difference between gas 
and radiation temperature, and the gas 
decouples earlier from radiation. All that 
tends to increase the 21 cm signal at high 
redshift. For redshifts z < 40, the situ- 
ation is more complicated. The altered 
evolution of the fluctuations and the de- 
lay in the formation of luminous objects 
decrease the amplitude of the power spec- 
trum in total, but also the onset of efficient 
coupling via the Wouthuysen-Field effect 
can be shifted significantly. 

In this context, one might wonder whether 
21 cm observations can actually probe ad- 
ditional small-scale power, which is pre- 



netic fields ( 


Kim, Olindo & Rosner 


1996 




Sethi & Subramanianll2005l: 


Tashiro & Sueivama 


2006). However, as discussed for ins 


tance 
), the 


by 


Furlanetto. Oh & Pierpaoli ( 


2006 



first 21 cm telescopes will focus on rather 
larger scales of the order 1 Mpc. The con- 
tribution to the power spectrum from pri- 
mordial magnetic fields depends on the 
assumed power spectrum for the mag- 
netic field. Regarding the formation of 
additional structures, the most significant 
case is a single-scale power spectrum. In 
this case, the contribution to the matter 
power spectrum c an be evaluated ana- 



lytically, yielding ( 


Kim. Olindo & Rosner 


1996 




Sethi & Subramanian 


2005) 



P(k) 



(30) 



for k < 2k max , where k max is given in 
Eq. [H For interesting field strengths of 
the order of 1 nG, we thus expect mod- 
ifications on scales of the order 1 kpc. 
For weaker fields, k^ ax is shifted to larger 
scales, but the power spectrum decreases 
with Bq. Detecting such changes in the 
matter power spectrum is thus certainly 
challenging. However, the calculation of 
the filtering mass indicates that the bary- 
onic gas will not collapse in minihalos of 
small masses when primordial fields are 
present. Star formation and HII regions 
from UV feedback are thus limited to more 
massive systems. This might be an addi- 
tional way to distinguish the case with and 
without magnetic fields. 

6. Discussion and outlook 

In the previous sections, we have pre- 
sented a semi-analytic model describing 
the post-reionization universe and reion- 
ization in the context of primordial mag- 
netic fields, and calculated the conse- 
quences for the mean 21 cm brightness 
fluctuation and the large-scale power spec- 
trum, formation of We identify two 
regimes in which primordial magnetic 
fields can influence effects measured with 
21 cm telescopes. At low redshifts, primor- 
dial magnetic fields tend to delay reioniza- 
tion and the build-up of a Lyman a back- 
ground, thus shifting the point where the 
signal is at its maximum, and changing the 
amplitude of the 21 cm power spectrum. 
The first 21 cm telescopes like LOFAR 
and others will focus mainly on the red- 
shift of reionization and can thus probe the 
epoch when a significant Lyman a back- 
ground builds up. As our understanding of 
the first stars increases due to advances in 
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theoretical modeling or due to better ob- 
servational constraints, this may allow us 
to determine whether primordial magnetic 
fields are needed to delay the build-up of 
Lyman a photons or not. As demonstrated 
above, comoving field strengths of the or- 
der 1 nG can delay the build-up of a Ly- 
man a background by Az ~ 10, which 
is significantly stronger than other mech- 
anisms that might delay the formation of 
luminous objects. Lyman Werner feed- 
back is essentially self-regulated and never 
leads to a significant suppress i on of star 
formation fl Johnson et all I2007L l2008h . X- 
ray feedback from miniquasars is strongly 
constrained from the observed soft X-ray 



backg round (ISalvaterra. Haardt fc Ferrara 
20051 ) . and significant heating from dark 



matter decay or annihilation would be 
accompanied by a significant amount of 
secondary ionization, resulting in a too 
large optical depth. On the other hand, 
an important issue is the question regard- 
in g the first sources of light. As shown 



in 



Schleicher. Banerjee fc Klessen 



massive Pop. Ill stars are needed to pro- 
vide the correct reionization optical depth. 
On the other hand, an additional popu- 
lation of less massive stars with an IMF 



according to (jScalol Il998l ; iKroupal 12002 
Chabrierl 120031 ) might be present. Such a 



population would emit more photons be- 
tween t he Lyman a and lin e per stellar 
baryon (ILeitherer et al.l Il999l ) , and could 
thus shift the build-up of a Lyman a back- 
ground to an earlier epoch. The onset of ef- 
ficient coupling via the Wouthuysen-Field 
effect thus translates into a combined con- 
straint on the stellar population and the 
strength of primordial magnetic fields. 



A further effect occurs at high redshifts, 
where the additional heat input from mag- 
netic fields due to ambipolar diffusion and 
the decay of MHD turbulence increases 
the 21 cm signal. As gas decouples ear- 
lier from the radiation field, the difference 
between gas and radiation temperature is 
larger and collisions are more effective in 
coupling the spin temperature to the gas 
temperature. The determination of the 21 
cm signal from this epoch is certainly chal- 
lenging, as the foreground emission corre- 
sponds to temperatures which are higher 
than the expected 21 cm brightness tem- 
perature by several orders of magnitude. 
However, as pointed out in other works 
flDi Matteo et all 120021; bh fe Mackl 120031: 



Zaldarriaea. Furlanetto fc H ernquist 



2004; 



ISethil 120051 : IShchekinov fc Vasilievl 120070 . 
the foreground emission is expected to be 
featureless in frequency, which may allow 
for a sufficiently accurate subtraction. In 
this context, it may also to help to fo- 
cus on the frequency gradient of the mean 
brightness temperature, rather than the 21 
cm brightness temperature itself. Upcom- 
ing long-wavelength experiments such as 
LOFAR, 21CMA (former PAST)|, MWA0, 
LWAI and SKA© may thus detect the addi- 
tional heat from primordial magnetic fields 
in the neutral gas, or otherwise set new 
upper limits on primordial magnetic fields, 
perhaps down to Bq ~ 0.1 nG. Like the 21 
cm transition of hydrogen, rotational and 
ro-vibrational transitions of primordial 
molecules may create i nteresting signa- 
tures in the CMB as well dSchleicher et al. 



3 http://wcb. phys.cmu.edu/ past 

4 http://web. haystack. mit.edu/arrays/MWA/index. html 
5 http: / /lwa. unm.edu 
6 http:/ /www. skatelescope.org 
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20081 ). which may provide a further test 
of the thermal evolution during the dark 
ages. 

We thank Benedetta Ciardi, Paul Clark, 
Christoph Federrath, Andrea Ferrara, Daniele 
Galli, Simon Glover, Thomas Greif, An- 
tonella Maselli and Francesco Palla for 
many interesting discussions on the topic. 
DRGS thanks the Heidelberg Graduate 
School of Fundamental Physics (HGSFP) 
and the LGFG for financial support. DRGS 
further thanks the MPA Garching for fi- 
nancial support during a research visit, 
and especially Benedetta Ciardi and An- 
tonella Maselli for being very kind host 
researchers. The HGSFP is funded by the 
Excellence Initiative of the German Gov- 
ernment (grant number GSC 129/1). RB is 
funded by the Emmy Noether grant (DFG) 
BA 3607/1. RSK thanks for support from 
the Emmy Noether grant KL 1358/1. All 
authors also acknowledge subsidies from 
the DFG SFB 439 Galaxies m the Early 
Universe. We thank the anonymous ref- 
eree for valuable comments that improved 
our manuscript. 

REFERENCES 

Abel, T., Bryan, G. L., Norman, M. L., 
2002, Science, 295, 93 

Allison, A. C, Dalgarno, A., 1969, AJ, 158, 
423 

Banerjee, R., Jedamzik, K., 2003, 
Phys. Rev. Lett., 91, 251301-1 

Banerjee, R., Jedamzik, K., 2004, 
Phys. Rev. D, 70, 123003 

Barkana, R., Loeb., A., 2001, Phys. Rep., 
349, 125 



Barkana, R., Loeb, A., 2005, MNRAS, 363, 
L36 

Barkana, R., Loeb, A., 2005, ApJ, 624, L65 

Barkana, R., Loeb, A., 2005, ApJ, 626, 1 

Becker, R. H. et al, 2001, AJ, 122, 2850 

Bharadwaj, S., Ali, S. S., 2004, MNRAS, 
352, 142 

Bland-Hawthorn, J., Maloney, P. R., 1999, 
ApJ, 510, L33 

Bland-Hawthorn, J., Maloney, P. R., 2001, 
ApJ, 550, L231 

Bromm, V., Kudritzki, R. P., Loeb, A., 
2001, ApJ, 552, 464 

Bromm, V., Larson, R., 2004, ARA&A, 42, 
79 

de Bruyn, A. G., Fender, R. P., Kuijpers, 
J. M. E., Miley, G. K., Ramachandran, 
R., Rttgering, H. J. A., Stappers, B. W., 
van de Weygaert, M. A. M., van Haar- 
lem, M. P., 2002, Exploring the Uni- 
verse with the Low Frequency Array. A 
scientific case. 

Chen X., Miralda-Escude J., 2004, ApJ, 
602, 1 

Chen, X., Miralda-Escude, J., 2006, eprint 
|arXiv: astro-ph/060543~9| 

Chabrier, G., 2003, PASP, 115, 763 

Choudhury, T. R., Ferrara, A., 2005, MN- 
RAS, 361, 577 

Christensson, M., Hindmarsh, M., Bran- 
denburg, A., 2001, Phys. Rev. E, 64, 
56405 



18 



Chuzhoy, L., Alvarez, M. A., Shapiro, P. 
R., 2006, ApJ, 648, LI 

Chuzhoy, L., Shapiro, P. R., 2006, ApJ, 
651, 1C 

Ciardi, B., Bianchi, S., Ferrara, A., 2002, 
MNRAS, 331, 463 

Ciardi B., Salvaterra R., 2007, MNRAS, 
381, 1137 

Clark, P. C, Glover, S. C. O., Klessen, R. 
S., 2008, ApJ, 672, 757 

Cowling, T. C, 1956, MNRAS, 116, 114 

Deharveng, J.-M., Buat, V., Le Brun, V., 
Millard, B., Kunth, D., Shull, J. M., 

Gry, C. A&A, 375, 805 

Di Matteo, T., Perna, R., Abel, T., Rees, 
M. J., 2002, ApJ, 564, 576 

Dove, J. B., Shull, J. M., Ferrara, A., 2000, 
ApJ, 531, 846 

Fernandez-Soto, A., Lanzetta, K. M., 
Chen, H.-W., 2003, MNRAS, 342, 1215 

Field, G. B., 1958, Proc. I.R.E., 46, 240 

Fujita, A., Martin, C. L., Mac Low, M.-M., 
Abel, T., 2003, ApJ, 599, 50 

Furlanetto, S. R., Oh, S. P., Briggs, F. H., 
2006, Physics Reports, 433, 181 

Furlanetto, S., Oh, S. P., Pierpaoli, E., 
2006, Phys. Rev. D, 74, 103502 

Furlanetto, S. R., Pritchard, J. R., 2006, 
MNRAS, 372, 1093 

Furlanetto, S. R., Furlanetto, M. R., 2007, 
MNRAS, 374, 547 



Giallongo, E., Cristiani, S., D'Odorico, S., 
Fontana, A., 2002, ApJ, 568, L9 

Gnedin, N. Y., Hui, L., 1998, MNRAS, 
296, 44 

Gnedin, N. Y., 2000, ApJ, 542, 535 

Goldman, S. P., 1989, Phys. Rev. A., 40, 
1185 

Greif, T. H., Johnson, J. L., Klessen, R. S., 
Bromm, V., 2008, MNRAS, in press 

Haiman, Z., Loeb, A., 1997, ApJ, 483, 21 

Heckman, T. M., Sembach, K. R., Meurer, 
G. R., Leitherer, C, Calzetti, D., Mar- 
tin, C. L., 2001, ApJ, 558, 56 

Hirata, C. M., 2006, MNRAS, 367, 259 

Hummer, D. G., Storey, P. J., 1987, MN- 
RAS, 224, 801 

Hurwitz, M., Jelinsky, P., Dixon, W. V. D., 

1997, ApJ, 481, L31 

Inoue, A. K., Iwata, I., Deharveng, J.-M., 
Buat, V., Burgarella, D., 2005, A&A, 
435, 471 

Jedamzik, K., Katalinic, V., Olinto, A. V., 

1998, Phys. Rev. D, 57, 3264 

Johnson, J. L., Greif, T. H., Bromm, V., 

2007, ApJ, 665, 85 

Johnson, J. L., Greif, T. H., Bromm, V., 

2008, MNRAS, 388, 26 

Kim, E.-J., Olinto, A. V., Rosner, R., 1996, 
ApJ, 468, 28 

Komatsu, E., et al., 2008, ApJS, submitted 



19 



Koyama, K., Petre, R., Gotthelf, E. V., 
Hwang, U. Matsuura, M., Ozaki, M., 
Holt, S. S., 1995, Nature, 378, 255 

Kroupa, P., 2002, Science, 295, 82 

Leitherer, C, et al, 1999, ApJS, 123, 3 

Loeb, A, Barkana, R., 2001, ARA&A, 39, 
19L 

Loeb, A., Zaldarriaga, M., 2004, 
Phys. Rev. Lett., 92, 211301 

Machida, M. N., Omukai, K., Matsumoto, 
T., Inutsuka, S.-i., 2006, ApJ, 647, LI 

Madau P., Meiksin A., Rees M. J., 1997, 
ApJ, 475, 429 

Malkan, M., Webb, W., Konopacky, Q., 
2003, ApJ, 598, 878 

Miiller, W., Biskamp, D., 2000, 
Phys. Rev. Lett., 84, 475 

Naoz, S., Barkana, R., 2005, MNRAS, 362, 
1047 

Nolta, M., et al., 2008, ApJS, submitted 

Oh, S. P., Haiman, Z., 2002, ApJ, 569, 558 

Oh, S. P., Mack, K. J., 2003, MNRAS, 346, 
871 

Olesen, P., 1997, Phys. Lett. B, 398, 321 

Omukai, K., Tsuribe, T., Schneider, R., & 
Ferrara, A. 2005, ApJ, 626, 627 

Omukai, K., Schneider, R., Haiman, Z., 
2008, submitted to ApJ 

Pequignot, D., Petitjean, P., Boisson, C, 
1991, A&A, 251, 680 



Pinto, C, Galli, D., 2008, A&A, submitted 

Press, W. H., Schechter, P., 1974, ApJ, 
187, 425 

Pritchard, J. R., Furlanetto, S. R., 2006, 
MNRAS, 367, 1057 

Purcell, E. M., Field, G. B., 1956, ApJ, 
124, 542 

Salvaterra, R., Haardt , F., Ferrara, A., 
2005, Monthly Notices of the Royal As- 
tronomical Society: Lett ers, 362, L50 

Sasaki, S., Takahara, F., 1993, PASJ, 45, 
655 

Scalo, J., 1998, in ASP Conf. Ser. 142, 
The Stellar Initial Mass Function, ed. 
G. Gilmore & D. Howell (San Francisco: 
ASP), 201 

Schleicher, D. R. G., Banerjee, R., Klessen, 
R. S., Phys. Rev. D, accepted (eprint 
larXiv:0807.3802ft 

Schleicher, D. R. G., Galli, D., Palla, F., 
Camenzind, M., Klessen, R., Bartel- 
mann, M., Glover, S. C. O., 2008, A&A, 
490, 521 

Schneider, R., Salvaterra, R., Ferrara, A., 
Ciardi, B., 2006, MNRAS, 369, 825 

Schneider, R., Salvaterra, R., Choudhury, 
T. R., Ferrara, A., Burigana, C, Popa, 
L. A., 2008, MNRAS, 384, 1525 

Seager, S., Sasselov, D., Scott, D., 1999, 
ApJ, 523, LI 

Seager, S., Sasselov, D., Scott, D., 2000, 
ApJ, 128, 407 



20 



Seshadri, T. R., Subramanian, K., 2001, 
Phys. Rev. Lett., 87, 101301 



Wong, W. Y., Moss, A., Scott, D, 2007, 
MNRAS, 386, 1023 



Sethi, S. K., Nath, B. B., Subramanian, 
K., 2008, MNRAS, accepted 



Wood, K., Loeb, A., 2000, ApJ, 545, 86 



Wouthuysen, S. A., 1952, AJ, 57, 31 



Sethi, S. K., Subramanian, K., 2005, MN- 
RAS, 356, 778 



Zaldarriaga, M., Furlanetto, S. R., Hern- 
quist, L., 2004, ApJ, 608, 622 



Sethi, S. K., 2005, MNRAS, 363, 818 



Zygelman, B., 2005, ApJ, 622, 1356 



Shapiro, P. R., Giroux, M. L., 1987, ApJ, 
321, L107 

Shchekinov, Yu. A., Vasiliev, E. O., 2007, 
MNRAS, 379, 1003S 

Shiromizu, T., 1998, Phys. Lett. B, 443, 
127 

Shull, J. M., van Steenberg, M. E., 1985, 
ApJ, 298, 268 

Silk, J., Langer, M., 2006, MNRAS, 371, 
444 

Smith, F. J., 1966, Planet. Space Sci., 14, 
929 

Steidel, C. C, Pettini, M., Adelberger, K. 
L., 2001, ApJ, 546, 665 

Subramanian, K., Barrow, J. D., 1998, 
Phys. Rev. D, 58, 83502 

Tashiro, H., Sugiyama, N., 2006, MNRAS, 
372, 1060 

Tashiro, H., Sugiyama, N., 2006, MNRAS, 
368, 965 

Whalen, D., Abel, T., Norman, M. L., 
2004, ApJ, 610, 14 

Wise, J. H., Cen., R., 2008, eprint 

arXiv:0808.2477 This 2-column preprint was prepared with the 



AAS IATgX macros v5.2. 



21 



